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1.  Introduction 


Infrasound  is  the  area  of  acoustics  that  is  below  the  audible  region  of  human  hearing,  starting  at 
15-20  Hz.  For  Army  infrasound  applications,  the  frequency  range  of  interest  is  in  the  1-20  Hz 
band.  Unlike  audible  acoustic  detection  ranges  typical  of  5-20  km  for  most  sources,  infrasound 
has  the  ability  to  propagate  extremely  long  ranges  of  100-500  km  and,  in  some  cases,  much 
further.  For  infrasound,  one  of  the  primary  sources  of  noise  is  wind  noise.  Wind  noise,  generated 
by  the  passage  of  turbulent  air  over  the  microphone,  is  often  suppressed  by  a  windscreen 
comprised  of  four  to  six  nominally  20-ft  lengths  of  porous  hose  radiating  out  from  the 
microphone.  This  windscreen  has  limited  the  ability  to  collect  infrasonic  signatures  in  covert 
areas  and  to  deploy  arrays  in  areas  where  there  is  limited  available  space.  For  both  military 
applications  and  environmental  remote  sensing  studies,  it  lacks  portability.  Physically,  the  porous 
hoses  are  acting  as  integrated  spatial  filter.  While  there  are  other  developmental  windscreens, 
none  have  been  fully  demonstrated  and  documented.  This  report  presents  the  results  of  a 
theoretical  and  experimental  examination  of  alternative  windscreens,  specifically  porous  and 
non-porous  polymer  fabrics  stretched  in  dome  configurations. 

For  audible  frequency  acoustics,  the  typical  windscreen  is  a  sphere  made  of  a  porous  material 
that  integrates  the  small-scale  turbulent  fluctuations,  leaving  a  near  laminar  flow  at  the  center  of 
the  windscreen.  From  Kolmogorov  (i),  the  turbulence  velocity  spectra  (figure  1)  are  separated 
into  frequency  ranges.  These  ranges  are  typically  expressed  as  three  spatial  scales  of  turbulence 
source  region  (large  scale),  inertial  subrange  (intermediate  scale),  and  dissipation  region  (small 
scale).  For  infrasound,  the  region  of  interest  for  wind  noise  generation  is  the  source  region,  and 
the  low  frequency  part  of  the  inertial  sub-range,  which  is  comprised  of  large  eddies  with  length 
scales  of  tens  of  meters  to  kilometers.  Unlike  higher  frequency  turbulence,  the  infrasound  region 
is  not  isotropic  and  not  homogeneous. 


Figure  1 .  The  Kolmogorov  turbulence  spectrum. 

In  the  early  2000s,  Liszka  (2),  Hedlin  and  Raspet  (3),  and  Bedard  et  al.  (4)  separately  showed 
that  wind  barriers  could  be  used  for  infrasonic  wind  noise  reduction.  Wind  barriers  break  up 
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turbulence  and  effectively  perform  a  spatial  average  on  turbulence  scales  smaller  than  the  size  of 
the  barrier.  In  these  studies,  a  large  wooden  structure,  such  as  in  figure  2a,  was  constructed 
around  the  infrasound  sensor  or,  in  Bedard’s  case,  a  porous  screen  fence  was  used  (figure  2b). 
The  results  of  these  studies  showed  equivalent  wind  noise  reduction  levels  as  the  ported  pipes 
being  used  at  the  time  for  Infrasound  Monitoring  Stations  (IMS).  The  research  has  continued  to 
study  different  designs  to  include  size,  materials,  and  extensions,  such  as  the  serrated  edges  first 
used  by  Bedard  to  break  up  the  turbulent  flow  coming  over  the  top  of  the  wind  barrier.  These 
ideas  and  designs  were  some  of  the  inspirations  behind  the  design  of  the  fabric  dome  used  in  this 
research  project. 


Figure  2.  Wind  barrier  designs  used  by  Liszka,  Hedlin,  and  Raspet  (a)  and  Bedard  (b). 

The  following  section  presents  a  theoretical  development  of  the  wind  noise  production  in  three 
regions  inside  and  outside  an  infrasound  windscreen.  Experimental  procedures  are  given  in 
section  3.  An  analysis  of  experimental  and  theoretical  results  is  given  in  section  4.  Section  5 
offers  some  concluding  remarks. 


2.  Theory 


The  one-dimensional  turbulent  kinetic  energy  spectrum  is  given  by  (5) 

E(k)  —  Cs2^3k~5^3  ,  (1) 

where  £  is  the  energy  dissipation  rate  of  turbulent  energy  per  unit  mass,  C  is  the  Kolmogorov 
constant  for  a  one-dimensional  spectrum,  which  typically  ranges  from  1.4  to  1.8,  and  k  is  the 
wavenumber.  The  energy  dissipation  rate  is  calculated  using 

£  =  +  Or  +  °"w)  ,  ( 2 ) 

where  a2,  and  a2  are  the  standard  deviations  of  the  it,  v,  and  w  components  of  the  wind 
field,  respectively.  Figure  3  shows  the  turbulent  kinetic  energy  spectrum  for  the  three  wind 
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conditions  used  in  the  present  analysis  of  the  wind  noise  measurements.  It  shows  that  as  the  wind 
speed  increases,  the  variance  in  the  wind  field  generally  increases,  providing  more  turbulent 
energy.  This  results  in  an  increase  in  the  wind  noise  detected  by  the  infrasound  microphone. 


Figure  3.  The  turbulent  kinetic  energy  spectrum  for  low, 
medium,  and  high  wind  conditions. 


For  low  frequency  acoustics,  the  wind  noise  contributions  due  to  turbulence  have  been  shown  to 
be  divided  into  turbulence-sensor,  turbulence-turbulence,  and  turbulence-mean  shear  interactions 
(<5).  Recently,  these  models  have  been  extended  to  include  contributions  in  the  source  region  of 
the  turbulence  spectrum  (7-9).  The  stagnation  pressure  is  a  result  of  the  pressure  imparted  on  an 
object  due  to  the  deflection  of  the  wind  around  the  object.  The  fluctuations  within  the  wind  field 
will  result  in  a  fluctuating  stagnation  pressure  on  the  sensor.  The  stagnation  pressure  power 
density  ( Ps )  due  to  the  turbulence-sensor  interaction  is 


Ps(k  i) 


1.44u2C  (  1.451C2 

[1  +  (M)2]5/6  +  A  [1  +  0.1129 (M)2]5/6  ' 


(3) 


where  ki  is  the  wavenumber  in  the  wind  direction  (k1  —  2nf  /u)  and  u  is  the  average  wind 
speed.  The  fitting  parameters  C  and  A  are  derived  from  the  one-dimensional  downwind  velocity 
spectrum  in  a  free  jet 


Vii(.ki)  [!  +  (/qA)2]5/6 


(4) 


using  ki  calculated  from  the  above  and  fitting  to  the  measured  velocity  spectrum.  This  can  be 
used  based  on  Taylor’s  frozen  turbulence  hypothesis,  which  states  the  turbulent  wind  field  is 
frozen  in  time  and  is  transported  past  the  sensor  at  the  mean  wind  speed,  u.  This  is  considered 
valid  since  the  lifetime  of  a  turbulent  eddy  is  long  compared  to  the  time  of  passage  by  the  sensor. 
Stagnation  pressure  is  the  pressure  measured  at  the  zero  velocity  position  on  a  rigid  sphere  and  is 
independent  of  the  radius.  This  formulation  developed  by  Raspet  et  al.  ( 8 )  extends  the  previous 
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formulation  for  the  stagnation  pressure,  which  was  only  valid  in  the  inertial  sub-range,  down  into 
the  source  region  of  the  turbulence  spectrum. 


The  next  source  of  noise  is  a  result  of  pressure  anomalies  generated  due  to  the  interaction 
between  turbulent  eddies.  This  effect  is  referred  to  as  the  turbulence-turbulence  interaction  and  is 
independent  of  the  presence  of  an  object.  For  the  inertial  sub-range,  this  noise  is  the  dominant 
source.  The  pressure  power  spectral  density  for  the  turbulence-turbulence  interaction  (Pt)  is 

C2  1 


Pt(k  r)=  0.811 


A  [1  +  0.1792(M)2]7/6  ' 


(5) 


This  form  of  the  equation  has  the  modifications  made  by  Raspet  et  al.  (8)  to  extend  it  into  the 
source  region,  which  is  the  dominant  source  of  turbulence  for  infrasound.  For  low  wavenumbers 
(large  scale  turbulence),  the  turbulence-turbulence  interaction  is  a  constant. 


Finally,  the  vertical  gradient  of  the  wind  field  near  the  surface  of  the  earth  dampens  turbulent 
eddies  being  generated  in  this  region  and  results  in  a  noise  known  as  the  turbulence-mean  shear 
interaction.  The  pressure  power  spectral  density  for  the  turbulence-mean  shear  interaction  (Pm)  is 


Pmik  i) 


7.38  CK2 


A2  (k^)5/3 

[1  +  1.622(M)2]8/3  ’ 


(6) 


where  K  is  the  gradient  of  the  average  wind  speed  at  the  observation  height.  Unlike  for  the  first 
two  interactions,  the  spectrum  at  low  wavenumbers  is  not  constant.  Instead,  the  spectrum 

5/3 

increases  as  kx  '  ,  peaks  just  before  the  transition  to  the  inertial  sub-range,  and  then  decays  as 


Based  on  work  by  Raspet  et  al.  (5),  Yu  et  al.  (7,  9,  10),  and  Abbott  et  al.  (11-13),  a  model  was 
developed  to  determine  the  pressure  at  the  center  of  a  wind  barrier  or  wind  screen.  The  pressure 
(Pc)  can  be  calculated  by  summing  the  pressure  contributions  due  to  atmospheric  turbulence  from 
each  of  the  three  regions  of  the  wind  screen:  inside  (Pi),  boundary  (P2),  and  outside  (P3).  These 
regions  are  shown  in  the  sketch  in  figure  4,  which  depicts  wind  flow  around  a  hemisphere  placed 
on  the  ground.  The  summed  pressure,  Pc,  is  expressed  as 

P2  =  P2  +  P2  +  P2  .  (7) 
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Figure  4.  Sketch  of  the  wind  flow  around  a  hemisphere  placed  on  the  ground. 


Starting  at  Region  1,  which  is  the  interior  of  the  porous  fabric  dome,  the  wind  field  will  be 
greatly  reduced  due  to  passage  through  the  porous  wind  barrier.  Abbott  (72)  suggests  modeling 
the  interior  wind  field  gradient  by  a  summation  of  exponentials 


(8) 


where  AL  and  /?j  are  fitting  parameters,  and  z  is  the  height  above  the  surface.  It  is  assumed  the 
dome  will  eliminate  the  stagnation  pressure,  thus  leaving  the  residual  wind  noise  contributions  to 
be  the  turbulence-turbulence  and  turbulence-mean  shear  interactions. 


Yu  (10)  derives  the  pressure  (p)  due  to  wind  noise  at  the  center  of  the  wind  barrier  by  starting 
with  the  source  equation  for  the  pressure  fluctuations  due  to  the  turbulence-mean  shear 
interaction 


V2p(r,  t)  = 


2  Po 


s(z)dVz 

dx 


(9) 


Here  r  =  (x,  y,  z)  is  the  coordinate  vector  where  x  corresponds  to  the  flow  direction,  y  is 
transverse  to  the  flow  (x,  y  are  in  the  horizontal  plane),  and  z  is  vertical.  Additionally,  p0  is  the 
density  of  air,  s(z)  =  dUx/dz  is  the  vertical  gradient  of  the  mean  wind  velocity  along  the 
direction  of  flow,  and  Vz  is  the  vertical  turbulent  velocity  component.  Applying  the  boundary 
condition  that  the  pressure  is  zero  at  the  surface  and  at  infinity  gives 


00  oo 

p(z,  k,  n>)  =  —  f  e~K\z'~z\T(z' ,  k,  cS)dz'  -\ - e~KZ  f  e~KZT(z',k,o))dz', 

2k  J  2k  J 


(10) 


where  T (z1 ,  k,  oj)  is  the  Fourier  transform  of  the  source  function  T (x,  t)  =  —  2p°s^dVzi  an(j 

OX 

k  —  kxex  +  ky  ey  is  the  wave  vector  in  the  plane  parallel  to  the  bounding  surface  (horizontal 
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plane).  Using  the  exponential  form  for  the  wind  velocity  profile  shown  in  equation  8  and  the 
Fourier  transform  of  the  velocity  correlation  tensor  applied  for  z  —  0  yields 


Ip(0,k)|2  = 


2po^x 


n 

i=n,j=n 


Z  s‘s'  /  ^  dkz  J  e  (K+Pi)z' sin (kzz')  dz'  J  e  (K+P/)z  sin(/ezz)  dz , 


(11) 


i—0j—0  —od  0 

where  E  (/r)  is  the  velocity  spectrum  function  given  by 

55C  (/cA)4 

ffb')  = _ v  J _ 

K  J  18  [1  +  (/rA)2]17/6  ‘ 
Substituting  this  into  equation  1 1  results  in 


(12) 


IPi(0,  kx)  |2 

i=n,j=m 


440 CA4plkl 


9n 


1 

.  JJ 


dkzdky 


1=0, J=0 


[1  +  (kA)2]17/6 


li  IL 

j  e-(/c+Pi)zsi n(kzz)dz  j  e_(K+^)z  sin (kzz')  dz’ 


(13) 


for  the  interior  of  the  fabric  dome,  where  h  is  the  height  of  the  wind  barrier,  p0  is  the  air  density, 
and  ^  and  fy  are  fitting  parameters.  The  wavenumber  is  k  —  ^Jk'x  +  kz  +  ky,  and  the 
wavenumber  for  the  horizontal-plane  projection  is  k  =  y/k2  -I-  ky  . 


The  next  region  is  the  undisturbed  flow  impinging  on  the  fabric  dome.  In  this  region,  the  wind 

XL  /  2  \ 

velocity  profile  is  assumed  to  be  a  typical  logarithmic  profile,  u  =  —  In  I  — )  where  u*  is  the 

'  kv  \Zq/ 


friction  velocity,  kv  is  the  Von  Karman  constant  (0.41),  and  z0  is  the  surface  roughness.  The 
friction  velocity  can  be  calculated  using 

u*  =  (itw)2  +  (vw)2,  (14) 


where  (itw)  and  (vw)  are  shearing  stresses  at  the  surface  for  the  momentum  fluxes  of  the  wind 
field.  Like  in  Region  1,  the  wind  noise  contribution  in  Region  3  can  be  calculated  from  the 
measured  turbulence  and  velocity  profile.  Going  through  the  same  process  as  before  except  using 
a  logarithmic  profile, 

it(z)  -  a  In  (— ),  (15) 

\z0/ 


instead  of  an  exponential  profile  and  using  the  trigonometric  function  difference  relations  results 
in 
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00 
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where  a  is  a  scaling  parameter  determined  by  the  wind  profile.  Substituting  equation  12  gives 
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(17) 


The  final  region  is  Region  2,  which  contains  the  surface  of  the  porous  fabric  dome.  Here  the 
stagnation  pressure  and  the  turbulence-mean  shear  interactions  are  the  dominant  sources  to  the 
wind  noise.  Currently,  there  is  only  an  empirical  model  to  estimate  the  contribution  to  the  wind 
noise  for  this  region  (12).  The  model  will  estimate  the  wind  noise  based  on  the  mean  wind 
velocity,  wind  gradient,  and  correlation  length  inside/outside  the  dome,  and  on  the  measured 
incident  turbulent  field.  Starting  with  the  Green’s  function  solution  to  the  Poisson  equation  for  an 
incompressible  fluid,  the  pressure  fluctuation  at  the  surface  of  the  dome  can  be  calculated  (14). 
The  solution  for  the  instantaneous  static  pressure  is 


p(r,t )  = 


Po  f  q(r',  t) 

4n  J  | f'  —  r\ 


(18) 


where  q(r' ,  t)  is  the  source  term  in  a  free  shear  flow.  Applying  the  Fourier  transform  to  the 
Green’s  function  and  converting  to  spherical  coordinates  yields 


p2(0,  k ) 


Po  ' 

47T„ 


Q(r,9,(p,k ) 
r 


r2sinddrd0d<p , 


(19) 


where  Q  (r,  9,  <p,  k )  is  the  Fourier  transform  of  the  source  term,  which  can  be  written  as 

Q(r,9,<p,k)=  iku(r,9,4>,k)S(9,4))  .  (20) 


The  term  S(6,  <fi)  represents  the  wind  velocity  gradient  across  the  dome’s  surface  and  can  be 
written  as 

Au 

S(6,<p)  =  —  ,  (21) 

o 

where  A u  is  the  measured  difference  in  the  wind  velocity  on  each  side  of  the  dome  surface,  and 
5  is  the  boundary  layer  thickness  shown  in  figure  5.  The  term  iku(r,  9,  (p,  k )  is  a  product  of  the 
wavenumber,  with  the  wind  velocity  turbulence  spectrum  impinging  on  the  surface  of  the  dome. 
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Figure  5.  The  boundary  layer  developed  by  the  wind  flow  over  the  dome. 


The  wind  velocity  spectrum  F(/c)  is  determined  by  performing  a  curve  fit  of  the  measured 
spectrum  to  the  von  Karman  spectrum  for  a  free  jet 


F(/c)  = 


C 

[1  +  (FA)2]5/6  ‘ 


(22) 


Substituting  equation  20  into  equation  19  yields 

p2(0,k)  —  —  — ^  f  u(r,  d,  <p,k)S {9,  tp)  r  sind dr dd dtp  .  (23) 

An  J 

Assuming  that  the  pressure  at  the  center  of  the  dome  is  dominated  by  the  thin  surface  layer 
{dr  ->  5),  equation  23  may  be  reduced  to  the  surface  integral 

n 

2n  2 

p2(0,k)  —  —  ^°  °  J  J  u{r0, 9,  <p,k)S  {6,  (p)  sind  dd  dtp  , 

o  o 

where  r0  is  the  radius  of  the  dome. 

Finally,  the  power  spectrum  for  the  pressure  is  calculated  by  taking  the  product  of  p2  with  its 
complex  conjugate  giving 


lp2((U)l2 


k2PorQ2(52 
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n  n 
2n  2  2n  2 

S*{9' ,  <p')S(d,  0)(n*(fo, (p1,  k)u(r0,  d,  tp,  k))sind'  sind  dd'  dtp'  dd  dtp  , 

oooo 


(25) 


where  {u*{Tq,  9',  <p',  k)u(r0, 9,  tp,  k))  is  defined  as  the  cross-correlation  of  the  wind  fluctuations. 
Studies  conducted  by  Priestley  ( 15 )  shows  that  the  cross-correlation  can  be  rewritten  as  the 
velocity  power  spectrum  times  an  exponential  decay  or 
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(u*(tq,9'  ,  <p',  k)u(r0, 9,  <p,  k))  =  u*(/c)w(/c)e“a(k)r°V(e'“0)2+  sin2  9  (0'“  *)2 


(26) 


where  the  attenuation  coefficient,  a(k),  is  related  to  the  inverse  of  the  correlation  length. 
Priestley  (15)  and  Shields  (16)  showed  that  the  attenuation  coefficient  could  be  expressed  in  the 
form  a(k)  —  fikck  where  q  and  q  are  fitting  parameters  dependent  on  the  porosity  of  the 
material.  Table  1  shows  the  values  for  q  and  q  determined  by  Abbott  (12)  for  a  number  of  wall 
porosities. 


Table  1.  The  values  for  the  fitting  parameters  q  and  q  for  different  fabric  porosities. 


Porosity 

q 

q 

30% 

1.0304 

0.7674 

40% 

1.4142 

0.6545 

55% 

1.6976 

0.6076 

65% 

1.2786 

0.7601 

80% 

1.5166 

0.6681 

Substituting  equation  26  into  equation  25  results  in 


|p2(0,fc)l2  =  l  if 


7i  n 
2n  2  2ir  2 
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(27) 


-a(k)r0y/(9r-9)2+  sin2  9  (_<p' -  <p)2 


sin9'sm9d9'  dtp' , 


where  u(k)  can  be  expressed  by  fitting  the  wind  velocity  spectrum  to  the  von  Karman  spectrum 
for  a  free  jet  (equation  22)  which  yields 


lp2«U)l2  =  j  1 1  jV(8',0')S(fl.0) 


n  n 
2tt  2  2tt  2 
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oooo 
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(28) 


e-a(fe)r0V ( 9'-9)2+  sin* 9  (0'-  <p)2  sin0'sin0d0'd0'd0d0. 


The  final  term  to  expand  is  the  product  of  the  velocity  gradients  at  the  surface  of  the  dome, 
S*(9' ,  <p')S(9,  <p).  The  estimate  of  the  velocity  gradient  is  given  in  equation  21.  As  mentioned 
previously,  A u  represent  the  difference  in  wind  velocity  between  outside  and  inside  of  the  dome 
or  A u  —  itout  —  itjn.  The  terms  itout  and  itin  are  considered  to  be  uniform  and  unidirectional.  As 
a  result,  A u  is  uniform  around  the  dome  for  a  turbulent  boundary  layer  thickness  8  defined  by 
u  —  0.99itout .  To  model  the  turbulent  boundary  layer  profile,  the  one-seventh  power-law 
velocity  profile  will  be  used  giving 
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u  fr\1^7 

lt0ut  3 

Solving  for  the  boundary  layer  thickness  results  in 


=  r(lW 


Substituting  this  into  equation  21  gives 


Au  (  u  \ 

5(0,0)  =  —  - 

r  \unuJ 


for  the  velocity  gradient  at  the  surface  of  the  dome  where  r  is  the  normal  distance  from  the 
surface,  itout  is  the  outside  measured  upwind  of  the  dome,  Ait  is  the  difference  in  wind  velocity 
between  outside  and  inside  of  the  dome,  and  u  is  the  wind  velocity  at  distance  r  from  the  surface 
of  the  dome.  To  determine  the  ratio  it/itout,  Abbott  (12)  used  the  cosine  Fourier  series 


it  v-1 

- —  A0  +  >  An  cos (n0)  , 

Unnt  2—J 


where  A{)  and  An  are  coefficients  determined  from  fitting  the  equation  to  measurements.  Finally, 
substituting  equation  32  into  equation  31  yields 


Au  V-1 

5(0,0)  =  —  A0  +  y  '  An  cos (n0) 


Substituting  this  result  into  equation  28  gives  the  final  form  for  the  power  spectrum  for  the 
pressure  contribution  to  the  wind  noise  in  Region  2 


2n  2  2n  2 

2„2„2r2/^2  r  r  r  r  2 


|p2(0,/c)|2  = 


16n2 
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Finally,  to  compute  the  total  contributions  to  the  wind  noise  at  the  center  of  the  dome,  equations 
28,  13,  and  17  can  be  substituted  into  equation  7  and  converted  from  wavenumber  space  to 

2n 

frequency  space  by  multiplying  by  — ,  where  un  is  the  mean  wind  velocity  in  the  region.  This 

un 

process  leads  to  the  following  expressions  for  the  pressures  in  each  region 
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and  the  total  at  the  center  of  the  dome 

lpc(o,/)l2  =  |Pi(0,/)|2+  |p2(0,/)|2+  |p3(0,/)|2. 


■dz' , 


(38) 


According  to  Abbott  (12),  the  comparisons  made  with  a  variety  of  porosities  under  different 
wind  conditions  indicated  that  for  low-  to  mid-level  porosities,  the  noise  contribution  due  to  the 
wind  barrier  is  the  most  dominant,  while  the  noise  contribution  due  to  the  flow  inside  the  fence  is 
negligible.  Conversely,  for  the  high  porosities,  the  contributions  of  the  pressure  fluctuations  on 
the  barrier’s  exterior  are  negligible,  while  the  inside  and  outside  contributions  are  the  most 
dominant.  The  comparisons  also  show  that  the  dominant  noise  sources  at  the  center  of  the  wind 
fence  for  low  wavenumbers  are  due  to  the  outside  flow  and  the  interactions  at  the  surface  of  the 
wind  fence,  while  noise  contributions  from  the  flow  region  inside  the  fence  are  negligible.  This 
is  consistent  for  all  porosities.  The  dominant  noise  sources  for  high  wavenumbers,  or 
wavenumbers  where  the  wavelengths  are  smaller  than  the  size  of  the  wind  fence  height,  are  due 
principally  to  the  interactions  at  the  surface  of  the  fence  for  low  porosities,  while  for  high 
porosities  the  dominant  noise  source  is  due  to  the  flow  interactions  inside  the  wind  fence.  Figure 
6  shows  Abbott’s  measured  pressure  spectra  inside  and  outside  of  his  wind  fence  compared  to 
the  theoretical  predictions.  Figure  6a  shows  the  contributions  to  the  pressure  spectra  from  each  of 
the  three  regions,  and  figure  6b  shows  how  the  summation  of  the  predicted  pressure  spectra 
compared  to  the  measured  spectra  inside  and  outside  the  wind  fence.  Abbott’s  model  shows  very 
good  agreement  to  the  measured  spectrum  for  this  case.  The  largest  contribution  to  the  overall 
wind  noise  is  shown  in  figure  6a  to  be  from  Regions  2  and  3.  This  case  was  chosen  due  to  the 
porosity  of  the  material  being  close  to  one  of  the  porosities  used  in  this  study. 
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Figure  6.  Abbott’s  measured  and  predicted  pressure  spectra  for  a  wind  fence  with  an  open  top  and  30%  wall 
porosity. 


3.  Experimental  Procedures 


To  investigate  the  ability  of  three  fabric  domes  to  successfully  mitigate  wind  noise  in  the 
infrasonic  range,  the  following  procedures  were  employed.  At  a  U.S.  Army  test  facility,  the  three 
domes  were  placed  close  to  one  another  and  at  an  angle  to  the  prevailing  winds  so  that  turbulent 
eddies  generated  by  the  passage  of  wind  over  one  of  the  domes  would  not  artificially  increase  the 
wind  noise  at  a  downstream  dome  (figure  7).  Hyperion  IFS-3000  infrasound  sensors  (the 
numbered  H’s  in  figure  7b)  and  Chaparral  Model  2  infrasound  sensors  (the  numbered  C’s  in 
figure  7b)  were  placed  at  the  center  of  each  dome.  Exterior  to  the  domes  were  two  IFS-3000 
sensors  and  one  Chaparral  Model  2.  The  Chaparral  and  one  of  the  IFS-3000  sensors  were 
equipped  with  a  20-ft  radius  porous  hose  rosette.  The  remaining  IFS-3000  was  equipped  with  a 
high-frequency  shroud  (HFS)  (figure  7c)  infrasound  windscreen  and  will  be  referred  to  as  the 
HFS-only  sensor  in  subsequent  discussions.  The  three  IFS-3000  sensors  in  the  domes  were  also 
equipped  with  HFS  windscreens.  All  signals  from  the  infrasound  sensors  were  recorded  using 
RefTek  130  data  recorders  sampling  at  1  kHz.  Only  the  internal  filters  in  the  RefTek  analog-to- 
digital  converters  were  used  to  low-pass  filter  the  infrasound  signals  at  the  Nyquist  frequency. 
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Figure  7.  Photograph  (a)  and  overhead  sketch  (b)  of  the  experimental  configuration, 
and  IFS-3000  with  high-frequency  shroud  (HFS)  (c). 

To  monitor  the  three-dimensional  air  motion  interior  and  exterior  to  the  domes,  R.M.  Young 
model  2000  ultrasonic  anemometers  (the  numbered  S’s  in  figure  7b)  were  used:  one  in  each 
dome  and  four  outside  the  domes.  The  ultrasonic  anemometer  measures  the  it,  v,  and  w 
components  of  the  wind  velocity  and  sonic  temperature  at  a  sample  rate  of  33  Hz.  Three  of  the 
exterior  anemometers  were  arranged  in  descending  height  (166,  123,  and  93  cm,  labeled  by  S4 
through  S6,  respectively,  in  figure  7b)  to  be  able  to  monitor  the  vertical  distribution  of  turbulence 
(72)  and  one,  S7,  was  placed  at  the  same  height  as  S5  to  monitor  the  horizontal  variation  of  the 
turbulence  flow  incident  on  the  domes.  An  Airmar  150WX,  A1  in  figure  7b,  at  a  height  of 
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159  cm,  was  used  to  collect  horizontal  wind  speed,  wind  direction,  temperature,  barometric 
pressure,  and  relative  humidity  at  a  sample  rate  of  1  Hz.  The  entire  experimental  setup  collected 
data  for  approximately  six  weeks  in  order  to  maximize  the  chances  of  wide  variations  in 
meteorological  conditions  and  to  collect  a  variety  of  impulsive  signatures. 

The  dome  shape  chosen  for  this  investigation  is  based  upon  a  commercially  available  2  m  radius 
tent,  which  was  used  as  the  non-porous  dome  and  is  referred  to  as  the  Nylon  dome.  Evidence  has 
been  presented  (13)  that  a  porosity  of  approximately  40%  works  well  for  infrasound  windscreen 
materials.  As  such,  two  porous  fabrics  were  chosen  to  verify  the  observations  made  in  Abbott  et 
al.  (13);  one  fabric  had  an  approximate  porosity  of  10%  (Sunbrella  Sling)  and  another  had  an 
approximate  porosity  of  35%  (Phifertex).  Both  porous  fabrics  are  vinyl-coated  polymer  strands, 
but  they  have  been  woven  in  such  a  fashion  as  to  be  porous.  Figures  8a  and  b  show  close-up 
photographs  of  each  fabric.  The  difference  in  porosity  is  clearly  visible  in  the  figure.  In  the 
following  discussion,  the  domes  constructed  from  each  fabric  will  be  referred  to  by  the  trade 
name  of  the  fabric. 


Figure  8.  Fabrics  used  in  the  construction  of  the  porous  domes;  Sunbrella  Sling  (a)  and  Phifertex  (b). 


4.  Results  and  Analysis 


In  order  to  determine  the  best  data  to  use  in  quantifying  the  wind  noise  reduction  characteristics 
of  each  of  the  domes,  the  wind  speed  and  direction  reported  by  the  Airmar  sensor  were  observed 
to  identify  times  when  the  wind  direction  was  ±45°  from  a  line  perpendicular  to  the  line  joining 
the  three  domes,  figure  7b.  The  wind  speeds  during  these  times  were  then  used  to  further 
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separate  the  data  into  time  segments,  where  the  wind  speed  fell  into  three  regions:  less  than 
3  m/s  (low),  between  3  and  6  m/s  (medium),  and  greater  than  6  m/s  (high).  The  power  spectral 
density  (PSD)  of  each  infrasound  sensor  was  then  computed  using  Welch’s  average  periodogram 
method  (17).  Figure  9  (low  wind  conditions),  figure  10  (medium  wind  conditions),  and  figure  11 
(high  wind  conditions)  depict  the  PSD  of  each  of  the  wind  noise  reduction  schemes  (a)  and  the 
difference  between  each  PSD  and  the  PSD  of  the  HFS-only  sensor  (b). 


Figure  9.  PSDs  (a)  of  the  Nylon  (blue),  Phifertex  (green),  Sunbrella  Sling  (red),  porous  hose  (cyan),  and  HFS 

(purple)  wind  noise  reduction  methods  in  the  low  wind  condition  (speed  less  than  3  m/s).  The  difference 
between  the  PSDs  of  the  Nylon  (blue),  Phifertex  (green),  Sunbrella  Sling  (red),  porous  hose  (cyan) 
reduction  methods  and  the  HFS  are  shown  in  b. 

It  is  readily  observed  in  figure  9a  (low  wind  conditions  with  wind  speeds  less  than  3  m/s),  that 
below  roughly  0.4  Hz,  the  Nylon  (blue),  Phifertex  (green),  and  Sunbrella  Sling  (red)  domes  tend 
to  reduce  the  wind  noise  relative  to  the  HFS-only  sensor  (purple)  by  a  similar  amount  to  the 
porous  hose  (cyan).  Between  0.4  and  4  Hz,  the  Phifertex  and  porous  hose  perform  roughly  the 
same,  while  the  Sunbrella  Sling  further  reduces  the  wind  noise  by  5  dB/Hz.  In  the  same 
frequency  range,  the  reduction  due  to  the  Nylon  dome  decreases  to  approximately  1.75  Hz. 
Between  4  and  10  Hz,  the  reductions  by  all  of  the  methods  except  the  HFS  are  within  5  dB/Hz  of 
each  other.  Above  10  Hz,  the  reduction  due  to  the  domes  remains  the  same,  while  the  reduction 
due  to  the  porous  hose  continues  to  increase.  These  differences  between  the  reduction  methods 
are  made  more  apparent  when  the  PSD  of  the  HFS-only  sensor  is  subtracted  from  the  PSDs  of 
the  other  methods,  as  shown  in  figure  9b. 

Figure  10a  shows  PSDs  of  each  of  the  reduction  methods  in  the  medium  wind  condition, 

3-6  m/s.  The  color  designations  are  the  same  as  in  figure  9.  As  expected,  all  of  the  PSDs  in 
figure  10a  have  increased  in  magnitude,  and  some  of  the  features  mentioned  in  the  discussion  of 
figure  9a  are  becoming  more  apparent.  Specifically,  the  increase  in  the  PSD  of  the  Nylon  tent 
centered  at  roughly  1.75  Hz,  and  the  similarity  between  the  reduction  due  to  the  porous  hose  and 
the  Sunbrella  Sling  dome  is  more  apparent.  In  this  medium  wind  condition,  however,  the 
Sunbrella  Sling  dome  maintains  its  noise  reduction  capabilities  at  frequencies  above  10  Hz, 
while  the  Nylon  and  Phifertex  domes  and  the  porous  hose  exhibit  a  decrease  in  their  noise 
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reduction  capabilities.  The  Phifertex  dome  begins  to  significantly  deviate  from  the  Sunbrella 
Sling  dome  at  roughly  6  Hz.  This  is  reinforced  by  the  differences  shown  in  figure  9b,  where  the 
Phifertex  and  Sunbrella  Sling  domes  have  PSDs  within  5  dB/Hz  from  0.1  to  6  Hz,  and  the 
Sunbrella  Sling  dome  and  the  porous  hose  array  remain  within  5  dB/Hz  of  each  other  from 
0.1  Hz  and  10  Hz. 


Figure  10.  PSDs  (a)  of  the  wind  noise  reduction  methods  for  the  medium  wind  condition  (speeds  from  3-6  m/s) 
and  differences  between  the  PSD  of  each  method  and  the  PSD  of  the  HFS  (b).  Color  designations  are 
the  same  as  in  figure  7. 

Figure  11a  shows  PSDs  of  each  of  the  wind  noise  reduction  schemes  in  the  high  wind  condition, 
speeds  greater  than  6  m/s.  Between  0.1  and  1  Hz,  all  of  the  reduction  schemes  are  within 
5  dB/Hz  (figure  1  lb).  Above  1  Hz,  the  HFS-only  sensor  and  the  Nylon  dome  significantly 
deviate  from  the  remaining  three  reduction  schemes.  The  Phifertex  dome,  until  8  Hz,  stays 
within  3  dB/Hz  of  the  Sunbrella  Sling  dome.  For  the  porous  hose,  this  frequency  range  extends 
to  roughly  10.5  Hz. 


Figure  11.  PSDs  (a)  of  the  wind  noise  reduction  methods  for  the  high  wind  condition  (speeds  greater  than  6  m/s) 
and  differences  between  the  PSD  of  each  method  and  the  PSD  of  the  HFS  (b).  Color  designations  are 
the  same  as  in  figure  7. 
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Figure  12  shows  the  turbulence  spectra  outside  (blue)  and  inside  the  Nylon  (green),  Phifertex 
(cyan),  and  Sunbrella  Sling  (purple)  domes  for  each  of  the  wind  conditions:  low  (a),  medium 
(b),  and  high  (c).  The  red  line  in  each  plot  is  the  theoretical  Kolmogorov  -5/3  turbulent  kinetic 
energy  spectrum  for  the  incident  turbulent  field.  In  figure  12,  as  might  be  expected,  it  is  observed 
that  the  incident  turbulence  spectrum  increases  as  the  wind  speed  increases.  This  is  also  true  of 
the  turbulence  spectrum  inside  the  Phifertex  dome,  although  the  dome  maintains  an  order  of 
magnitude  reduction  in  the  turbulence  relative  to  the  incident  spectrum.  The  turbulence  spectra 
inside  both  the  Nylon  and  Sunbrella  Sling  domes  do  not  behave  in  as  predictable  a  fashion  as  the 
Phifertex  dome.  If  figures  12a  through  12c  are  compared,  it  is  noted  that  for  frequencies  above 
0.5  Hz  the  spectrum  inside  the  Sunbrella  Sling  dome  remains  relatively  constant  regardless  of 
exterior  wind  speed.  At  lower  frequencies,  the  turbulence  spectrum  increases  with  wind  speed. 
Throughout  figure  12,  the  turbulence  spectrum  inside  the  Sunbrella  Sling  dome  has  a  curved 
shape  the  cause  of  which  is  not  known  at  this  time.  In  figure  11a,  the  turbulence  spectrum  of  the 
Nylon  dome  linearly  decreases  at  a  slope  slightly  less  negative  than  the  theoretical  -5/3.  As  the 
wind  speed  increases,  figure  lib  and  11c,  the  turbulence  spectrum  inside  the  Nylon  dome 
assumes  a  curved  shape  similar  to  that  inside  the  Sunbrella  Sling  dome.  At  the  highest  wind 
speeds,  the  reason  for  the  odd  shape  of  the  lower  frequency  portion  of  the  spectrum  is  not 
known. 
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Figure  12.  Turbulence  spectra  outside  (blue)  and  inside  the  Nylon  (green),  Phifertex  (cyan),  and  Sunbrella 
sling  (purple)  domes  for  the  low-  (a),  medium-  (b),  and  high-wind  (c)  conditions.  The  red  line  in 
each  plot  is  the  theoretical  incident  turbulent  kinetic  energy  spectrum. 

Figure  13  shows  a  single  explosion  event  (a)  captured  by  all  five  infrasound  sensors,  and  the 
average  transfer  functions  (b)  of  each  wind  noise  reduction  method  compared  to  the  infrasound 
sensor  with  only  the  HFS  installed.  In  figure  13a,  the  waveforms  have  been  time-aligned  to 
better  show  the  differences  between  the  recordings.  It  is  apparent  in  figure  13a  that  the 
waveforms  recorded  under  the  HFS  (taken  to  be  the  reference)  and  the  Phifertex  and  Nylon 
domes  are  all  similar  although  the  recording  in  the  Nylon  dome  has  higher  peak  amplitude  than 
the  HFS.  Also  noted  is  a  small  decrease  in  amplitude  of  the  waveform  due  to  the  Sunbrella  Sling 
dome  and  a  large  decrease  in  amplitude  due  to  the  porous  hose  that  is  accompanied  by  an 
apparent  low-pass  filtering  of  the  waveform.  This  significant  decrease  in  amplitude  due  to  the 
porous  hose  potentially  results  in  a  decrease  in  detection  range  when  using  porous  hose  as  a  wind 
noise  reduction  method.  The  transfer  functions  in  figure  13b  were  calculated  by  taking  the 
transfer  functions  between  each  wind  noise  reduction  method  and  the  HFS -only  sensor  due  to 
four  similar  waveforms  and  averaging  the  results.  The  error  bars  in  figure  13b  are  at  approximate 
third-octave  center  frequencies.  From  5  Hz  to  approximately  75  Hz,  the  Phifertex  dome  agrees 
very  well  with  the  HFS-only  sensor,  while  the  Sunbrella  Sling  dome  attenuates  the  signal  by 
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approximately  1  dB.  Thus,  both  porous  domes  perform  exceptionally  well  at  preserving  the 
properties  of  a  signal  passing  through  the  dome.  The  transfer  function  of  the  Nylon  dome  below 
20  Hz  is  flat  and  agrees  well  with  the  HFS-only  sensor.  However,  above  20  Hz,  the  transfer 
function  of  the  Nylon  dome  is  greater  than  one,  implying  an  increase  in  amplitude  of  the  signal 
reaching  the  sensor  at  the  center  of  the  dome.  The  transfer  function  of  the  sensor  with  the  porous 
hose  array  is  roughly  3  dB  below  the  HFS-only  sensor  at  5  Hz  and  steadily  decreases  to  4  dB 
below  at  20  Hz.  Above  20  Hz,  the  response  falls  rapidly  to  approximately  12  dB  below  the  HFS- 
only  sensor  at  about  60  Hz.  This  shows  the  low-pass  filter  effect  of  the  porous  hose  alluded  to  in 
the  discussion  of  figure  13a. 


Figure  13.  Waveforms  of  a  small  explosion  (a)  and  the  average  transfer  functions  (b)  of  each  wind  noise  reduction 
method  relative  to  the  sensor  with  only  the  HFS  installed. 

Based  upon  the  wind  noise  spectra,  the  turbulence  spectra,  and  the  transfer  functions,  the  porous 
domes  have  demonstrated  a  distinct  ability  to  reduce  wind  noise  at  an  infrasound  sensor  and  to 
preserve  the  characteristics  of  a  signal  passing  through  the  wall  of  the  domes.  Of  the  two  porous 
domes,  the  Sunbrella  Sling  appears  to  have  the  best  overall  performance:  it  maximizes  the  wind 
noise  reduction  while  losing  only  a  small  amount  of  signal.  The  Phifertex  dome  clearly 
demonstrated  the  best  signal  preservation,  but  its  ability  to  reduce  the  wind  noise  reaching  the 
sensor  was  not  sufficiently  high  at  the  upper  end  of  the  typical  frequency  range  used  in  Army 
applications,  1-20  Hz.  The  non-porous  dome  introduced  too  many  artifacts,  as  discussed 
previously,  into  the  signals  to  consider  it  as  a  viable  alternative  to  porous  hose. 
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5.  Concluding  Remarks 


Two  porous  fabric  and  one  non-porous  fabric  domes  have  been  investigated  to  determine  their 
potential  as  alternative  methods  for  infrasound  wind  noise  reduction.  Transfer  functions  between 
a  reference  sensor  and  signals  measured  at  the  center  of  the  domes  or  at  the  center  of  a  porous 
hose  array  demonstrated  the  ability  of  the  porous  domes  to  preserve  the  information  content  in  a 
signal.  Further,  wind  noise  spectra  measured  at  the  center  of  the  domes,  when  compared  to 
spectra  measured  by  the  reference  sensor,  demonstrated  the  ability  of  the  domes  to  reduce  the 
wind  noise  to  a  level  comparable  to  that  of  a  porous  hose  array.  Of  the  two  porous  domes 
investigated,  the  dome  constructed  from  Sunbrella  Sling  fabric  exhibited  the  best  performance 
and  can  be  considered  as  a  recommended  alternative  to  porous  hose  rosettes. 
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